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f*^) , Abstract. Deep-inelastic e + p scattering data, taken with the HI detector at HERA, are used 

to investigate jet production over a range of four-momentum transfers 150 < Q 2 < 15000GeV 2 
Cn| I and transverse jet energies 5 < Ej < 50GeV. The analysis is based on data corresponding to an 

integrated luminosity of Jz?j nt = 65.4pb~' taken in the years 1999-2000 at a centre-of-mass energy 
y/s 319GeV. Jets are defined by the inclusive k, algorithm in the Breit frame of reference. 
Dijet and trijet jet cross sections are measured with respect to the exchanged boson virtuality 
and in addition the ratio of the trijet to the dijet cross section R^n is investigated. The results 
| are compared to the predictions of perturbative QCD calculations in next-to-leading order in the 

strong coupling constant a s . The value of a s (mz) determined from the study of R^/ 2 is 0C s (mz) = 

0.1175 ±0.0017(stat.)±0.0050(syst.)+°;°^(theo.). 

On 

O INTRODUCTION AND OBSERVABLES 

O ' 

. Deep-inelastic scattering (DIS) at HERA is a precision tool for studies of Quantum 

Chromo Dynamics (QCD). Jet cross sections at high transverse momenta are in particu- 
lar attractive, because perturbative calculations (pQCD) are precise and non-perturbative 
hadronisation effects are weak. In the past inclusive jet cross sections QjJ, |2|] , as well as 



Qh! multijet cross sections 10, HI] have been studied at HERA. 

The aim of the present analysis is to check quantitatively pQCD predictions for di- 
and trijet cross sections and to determine the value of the strong coupling at different 
values of the relevant hard scale. 

In the following, jets are defined by the inclusive k t cluster algorithm Hi in the Breit 
frame of reference. Accepted jets are required to have a transverse energy of more than 
5 GeV. These jets, boosted back to the laboratory frame, have to fulfill the pseudorapidity 
requirement — 1 < 7]i a b < 2.5 to ensure good detector acceptance. Events with at least 
two (three) accepted jets are assigned to the dijet (trijet) event sample. In order to ensure 
stability of the perturbative calculations, cuts on the invariant jet masses are applied: 
Mdijet > 25 GeV and M tr ij et > 25 GeV for the di- and trijet sample, respectively. 

DATA SAMPLE AND ANALYSIS METHODS 

The data this analysis is based on were taken with the HI detector in the years 1999- 
2000 and correspond to an integrated luminosity of J2?i nt = 65.4pb _1 . A DIS selection 
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FIGURE 1. NC dijet and trijet differential cross-sections, with respect to Q 2 , shown with NLO pQCD 
predictions including hadronisation corrections. The shaded bands show the effect of varying the renor- 
malisation/factorisation scale by a factor of two. 



in the phase space of four-momentum transfers 150 < Q 2 < 15000 GeV 2 and inelasticity 
0.2 < y < 0.6 leaves 5460 dijet- and 1757 trijet events. 

To account for limited detector acceptance and resolution, correction factors were 
applied, which were determined with the Monte Carlo event generators DJANGOH 0] 
(including ARIADNE |8]) and RAPGAP The cross sections were corrected for QED 
radiative effects with HERACLES ifToll . 

Relevant experimental uncertainties include: the electromagnetic energy scale, the 
hadronic energy scale, the scattering angle of the positron, a model uncertainty in the 
acceptance correction and the luminosity measurement, with the hadronic energy scale 
uncertainty being the dominant contribution. 



RESULTS 

Fig. [l] shows the differential di- and trijet cross sections as a function of Q 2 . The data 
are compared to a calculation carried out with NLOJET++ (TH| . using parton density 
functions of the proton from the CTEQ5M1 lfl2ll set and assuming a value of the strong 
coupling a s (mz) =0.118. The calculation includes matrix elements at next-to-leading 
order (NLO) for two and three partons in the final state, corresponding to orders of the 
strong coupling up to ff(a^). Hadronisation corrections are applied to the parton level 
prediction, determined with the help of DJANGOH and RAPGAP event samples. An 
uncertainty estimate for the prediction is provided by variations of the renormalisation 
and factorisation scale by a factor of two. The cross sections span several orders of 
magnitude, where a good description of the data by the prediction is observed. At the 
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FIGURE 2. Measured values of i 2 against Q 2 compared with a NLO pQCD prediction with hadroni- 
sation corrections. The dark shaded band shows the uncertainty associated with the hadronisation correc- 
tions, while the light shaded band shows the effect of varying the renormalisation/factorisation scale by a 
factor of two. 



highest Q 2 bin electroweak effects, which are not present in the calculation, cannot be 
neglected. Hence this point will not be used in the following a s fit. 

The ratio of the tri- and dijet cross section, R^n = Otrijet/°dijet is shown in Fig. EJ 
where again a good description by the perturbative calculation is observed. Based on this 
measurement, a fit of the strong coupling is performed: for each Q 2 value at which i? 3 / 2 is 
measured (but the highest one) the pQCD calculation is repeated with the five variations 
of the proton p.d.f.s, available within the CTEQ4A II 1 3fl set. This yields five predictions 
of i? 3/ / 2 as a function of a s (mz), to which a function of the form C\(X s (mz) +C20cj(mz) 
is fitted. Consequently, the function is used to obtain from the measured value of i? 3 / 2 
and its uncertainty the corresponding value and uncertainty of a s (mz). 

The fit results are shown as points in Fig. |3] From this diagram the running of the 
renormalised strong coupling is clearly evident. In addition an average value is obtained 
from this values by % 2 minimisation, shown as a diamond at the reference of Q 2 = mk . 
The presented result is well compatible with the world average from the PDG H14I1 . 
which is shown as a band for comparison. 



CONCLUSION 

Differential di- and trijet cross sections measured at high Q 2 have been shown. The 
distributions are well described by NLO pQCD with hadronisation corrections, except 
for the highest Q 2 data point, where electroweak corrections cannot be neglected. A fit 
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FIGURE 3. a s (mz) values from each Q 2 bin evolved to values at their respective values of Q 2 (points) 
using the two-loop solution of the renormalisation group equation. The averaged value of a s (mz), found 
using a % 2 minimisation fit, is shown at the far right (empty diamond). The inner error bars show the 
statistical errors alone and the outer error bars denote the quadratic sum of statistical and systematical 
errors. The evolution of the world average value of a s (mz) is shown as a shaded band. 



of the strong coupling yields 

a s (m z ) = 0.1175 ± 0.0017(stat.) ± 0.0050(syst)^;^(theo.) , 

well in agreement with the world average. Future improvements in the understanding of 
the hadronic energy scale will significantly reduce the experimental uncertainty. 
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